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ABSTRACT 

Observations with the Spitzer Space Telescope have recently revealed a significant population of 
high-redshift (z ~ 2) dust-obscured galaxies (DOGs) with large (rest-frame) mid-infrared to ultra- 
violet luminosity ratios. Due to their optical faintness, these galaxies have been previously missed 
in traditional optical studies of the distant universe. We present a simple method for selecting this 
high-redshift population based solely on the ratio of the observed mid-infrared 24/im to optical i?-band 
flux density. We apply this method to observations of the ~ 8.6 deg 2 Bootes Field of the NO AO Deep 
Wide-Field Survey, and uncover «2,600 DOG candidates (i.e., a surface density of 0.089 arcmin -2 ) 
with 24/im flux densities F 2 ^m > 0.3mJy and (R - [24]) > 14 (i.e., F v (2^m) / F V (R) > 1000). These 
galaxies have no counterparts in the local universe. They become a larger fraction of the population 
at fainter 24^m flux densities, increasing from 7±0.6% of sources at i*24/*m > 1 mJy to w 13 ± 1% 
of the population at ss 0.3 mJy. These galaxies exhibit evidence of both star-formation and AGN 
activity, with the brighter 24/xm sources being more AGN-dominated. Their mid-infrared spectral 
energy distributions range from power-laws (likely AGN-dominated at mid-IR wavelengths) to sys- 
tems showing a "bump", the latter likely resulting from the redshifted 1.6/im peak characteristic of 
most stellar populations. Using primarily the W. M. Keck Observatory and Spitzer, we have obtained 
spectroscopic redshifts for 86 objects within this sample, and find a broad redshift distribution which 
can be modeled as a Gaussian centered at z ~ 1.99 ± 0.05 and a(z) « 0.45 ± 0.05. The space density 
of this population is Sdog (i*24/«n > 0.3 mJy) = (2.82 ± 0.05) x 10~ 5 /iy Mpc -3 , similar to that of 
bright sub-millimeter-selected or UV-selected galaxies at comparable redshifts. These redshifts also 
imply very large luminosities, with a sample median vL u [8^,m) «4x lO 11 ^©, implying 8[im — 1mm 
luminosities of Lir^ 10 12_14 Lq for the population. The infrared luminosity density contributed by 
this relatively rare DOG population is log(Li R ) » 8.23±g;^g. This is « 60tfg% of that contributed by 
z ~ 2 ultraluminous infrared galaxies (ULIRGs, with Ltr > 10 12 L©), and suggests that our simple 
selection criterion effectively identifies a significant fraction of z ~ 2 ULIRGs. This IRLD is also 
s=s 26 ± 14% of the total contributed by all z ~ 2 galaxies, and comparable to that contributed by 
the luminous UV-bright star- forming galaxy populations at z ~ 2. We suggest that these DOGs are 
the progenitors of luminous (~ AL*) present-day galaxies and are undergoing an extremely luminous, 
short-lived phase of both bulge and black hole growth. They may represent a brief evolutionary phase 
between sub-millimeter-selected galaxies and less obscured quasars or galaxies. 
Subject headings: galaxies:formation-galaxies:high-redshift-galaxies:starburst 



1. INTRODUCTION 

Infrared-luminous sources (with Ltr > 10 L©) domi- 
nate the bright end of the bolometric luminosity function 
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of galaxies in the local universe (Soifer et al. 1986). Their 
role in the story of galaxy evolution is not yet understood, 
but it is likely that their prodigious luminosities are the 
product of an extremely active phase during which these 
systems form stars and / or grow their central black holes 
at a rapid rate (Sanders et al. 1988). Did all large galax- 
ies exhibit such a phase in their early formative stages, 
or is this active phase characteristic of only some of the 
most massive systems? To address these questions, we 
need to find and study the primary high-redshift galaxy 
populations contributing to the mid-IR and far-IR counts 
and understand their relation to the local galaxy popu- 
lation. 

We know a great deal about the local population of 
IR-luminous galaxies, due in large part to the pioneer- 
ing observations made with the Infrared Astronomical 
Satellite (IRAS; e.g., Soifer et al. 1987; Sanders et al. 
1988; Sanders & Mirabel 1996). Their mid-IR number 
counts show strong evolution (e.g., Elbaz et al. 1999; 
Chary & Elbaz 2001; Papovich et al. 2004; Chary et al. 
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2004; Marleau et al. 2004). Although rare locally, in- 
frared luminous galaxies become an increasing fraction 
of the galaxy population at high redshift, dominating the 
IR energy density at redshifts z ~ 1 (e.g., Franceschini 
et al. 2001; Lagache et al. 2004; Gruppioni et al. 2005; 
Le Floc'h et al. 2005; Perez-Gonzalez et al. 2005). 

Far less is known about their higher redshift counter- 
parts despite several ground-breaking studies. Prior to 
the launch of the Spitzer Space Telescope, our knowledge 
of high-z IR-luminous galaxies was derived from studies 
of a handful of the most extreme IRAS sources, con- 
straints on the redshift evolution out to z ~ 1 (largely 
due to ISO; e.g., Franceschini et al. 2001; Chary & Elbaz 
2001; Lagache et al. 2003), small samples of galaxies dis- 
covered at sub-millimeter wavelengths (e.g., Smail et al. 
2002; Blain et al. 2002; Chapman et al. 2004, 2005), and 
several galaxies discovered as a result of their extreme 
optical-to-near-infrared colors (the dusty 'Extremely Red 
Objects', or EROs, with R — K > 6; e.g., Hu & Ridgway 
1994; Dey et al. 1999; McCarthy 2004). 

The launch of the Spitzer Space Telescope has revolu- 
tionized studies of the dusty galaxy population. It has 
provided the first mid-IR spectroscopy of dust emission 
in high-redshift star-forming galaxies and AGN, and re- 
vealed populations of galaxies that are extremely faint at 
(observed) visible wavelengths (e.g., Houck et al. 2005; 
Yan et al. 2005, 2007, Higdon et al. and Desai et al., 
both in preparation). Many of these galaxies lie at red- 
shifts z > 1, and their mid-IR spectra reveal either strong 
silicate absorption features, suggestive of an AGN, or 
emission from polycyclic aromatic hydrocarbons (PAHs) 
which generally arise in photodissociation regions asso- 
ciated with star-forming molecular clouds. A significant 
fraction of these active high-redshift galaxies are under- 
luminous at rest-frame UV wavelengths, suggesting that 
the prodigious star-formation and / or AGN accretion 
luminosity is hidden by dust that reprocesses and radi- 
ates the bulk of the energy at far-IR wavelengths (e.g., 
Sanders & Mirabel 1996). 

In this paper, we present a simple and economical 
method for selecting these high-redshift dusty popula- 
tions, based solely on the ratio of the observed 24^,m to 
optical flux ratio. This selection, applied to a wide area 
survey, results in a large, significant population of ex- 
tremely luminous high-redshift galaxies for which there 
are no local counterparts. These are systems with un- 
usually red spectral energy distributions, which we inter- 
pret as being the result of extinction by large columns of 
dust. These galaxies appear to be largely absent from 
(rest-frame) UV-selected high-redshift galaxy samples. 
We have measured a large number of spectroscopic red- 
shifts for this population, and confirm that they lie at 
redshifts z ~ 2. Given their large mid-infrared lumi- 
nosities and implied far-infrared luminosities, these dust 
obscured galaxies may be undergoing an intense burst of 
star formation, accretion activity (onto a nuclear black 
hole) or plausibly both. We speculate that this popu- 
lation of objects, which may be selected simply on the 
basis of optical to mid-infrared flux ratio, are good can- 
didates for systems in the throes of formation of their 
stellar spheroids and nuclear black holes. Despite their 
rarity, this population contributes a significant fraction 
of the infrared luminosity density at z ~ 2. 



This paper is structured as follows. In §2, we briefly 
describe the observations from which the sample is drawn 
and discuss the selection criteria we use to isolate the can- 
didate high-redshift dust-obscured galaxy (hereinafter 
DOG) population. In §3, we present the main re- 
sults: the optical-to-mid-infrared color distribution of 
the DOGs, their spectral energy distributions, redshift 
distribution and rest-frame UV and IR luminosities. In 
§4, we discuss the nature of the DOG population, derive 
their space density and their contribution to the z m 2 
infrared luminosity density, and discuss their possible re- 
lationship to other z f=s 2 galaxy populations. We sum- 
marize our findings in §5. 

Throughout this paper we adopt a cosmology with 
H = 70 km s" 1 Mpc^ 1 , n m = 0.3, fi A = 0.7. Mag- 
nitudes quoted in this paper arc in Vega units, except 
where explicitly defined otherwise (e.g., in the discussion 
of rest-frame 2200A luminosities). 

2. OBSERVATIONAL DATA AND SAMPLE SELECTION 

2.1. Imaging Data 

The study described in this paper focusses on a sam- 
ple drawn from the mid-infrared and optical observa- 
tions of the Bootes Field of the NO AO Deep Wide-Field 
Survey (NDWFS 13 ; Jannuzi & Dey 1999, Jannuzi et 
al., in prep., Dey et al., in prep.). The NDWFS is a 
deep, ground-based, optical and near-infrared imaging 
survey of two 9.3 square degree fields, one in Bootes 
and one in Cetus, conducted using the 4m and 2.1m 
telescopes of the National Optical Astronomy Observa- 
tory. The survey reaches median 5<t point-source depths 
in the B W ,R,I and K bands of w27.1, 26.1, 25.4 and 
19.0 mag (Vega) respectively (see the NOAO Science 
Archive at http : //http : //archive . noao . edu/nsa/ for 
details). The data products for the NDWFS will be de- 
scribed elsewhere (Dey et al., in prep.; Jannuzi et al., 
in prep). The NDWFS astrometry is tied to the refer- 
ence frame defined by stars from the USNO A- 2 catalog 
(Monet et al. 1998). 

Spitzer has mapped 8.61 deg 2 of the NDWFS Bootes 
field at 24, 70 and 160/im using the Multiband Imaging 
Photometer for Spitzer (MIPS; Rieke et al. 2004) to la 
rms depths of 0.051mJy, 5mJy and 18mJy respectively. 
The data were reduced by the MIPS GTO team. De- 
tails of the survey, including the mapping strategy, data 
reduction and the resulting catalog, will be discussed 
elsewhere. The entire Bootes field was also mapped at 
3.6, 4.5, 5.8 and 8.0/xm using the Infrared Array Cam- 
era (IRAC; Fazio et al. 2004) on Spitzer, details of the 
IRAC observations may be found in Eisenhardt et al. 
(2004). Both the MIPS and IRAC data are astrometri- 
cally calibrated to the 2MASS astrometric frame. There 
are roughly 22,000 24/im sources in the Bootes field down 
to the 80% completeness limit of 0.3 mJy. 

The 24/zm source catalogs were matched to the opti- 
cal catalogs (after accounting for a small but significant 
relative offset of Aa = 0'.'38 and AS = 0'.'14 between 
the two frames). Optical photometry was determined 
for each 24/zm source using the images of the third data 
release (DR3) of the NDWFS. In order to measure ac- 
curate optical colors in small aperture photometry, we 
smoothed the optical NDWFS images so the delivered 

13 See also http://www.noao.edu/noaodeep/ 
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point spread function across the entire Bootes field was 
a 1.35" FWHM Moffat profile with a /3 parameter of 2.5. 
Using our own code, we measured 4" diameter aperture 
photometry at the 24 /j,m position of each source. The 
i?-band magnitudes of the brighter sources, i.e. with 
R < 25, are the 'AUTO' magnitudes derived by Source- 
Extractor (Bertin & Arnouts 1996). For fainter sources 
and for mid-IR sources without optical counterparts in 
the optical catalogs, photometric limits were measured 
directly from the optical stacked images in 4" diame- 
ter apertures. In cases where the 4" diameter aperture 
magnitude (or limit) was indeterminate or greater than 
the 50% completeness limit of our optical data, we re- 
placed the estimated magnitude with the brighter of the 
locally measured 3cr sky noise or the 50% completeness 
limit. The background and uncertainties of the photom- 
etry were determined by measuring the sigma-clipped 
mean and RMS of fluxes measured in sixteen 4" aper- 
tures located in a 45" radius circle around each object. 

2.2. Sample Selection 

Our sample is selected from the 8.140 deg 2 region 
of the Bootes field for which there exists both good 
Spitzer/MiPS 24/xm and KPNO i?-band imaging. Re- 
gions around bright, saturated stars were masked out. 
(R — [24]) colors (in the Vega system) 14 computed for 
all 24/im sources are shown as a function of 24/im flux 
density (i*24/im) m Figure 1. The points appear to 
spread toward redder colors at fainter 



F v l24nmJ (mJy) 



24/tm ■ 



In or- 



der to quantify this effect, we define a fiducial color of 
(R - [24]) = 14 corresponding to F^ m /F R > 982, cho- 
sen because this is redder than the color of most ultra- 
luminous infrared galaxies (ULIRGs; > 10 12 L Q ) at 
all redshifts (Figure lb). We measure the (differential) 
fraction of the source population redder than this limit as 
a function of -F24/im and find that the fraction increases 
to fainter flux densities, rising from 7±0.6% of the pop- 
ulation for all sources with i^/im > lmJy to 13±1% at 
Funm « 0.3mJy (Figure 2). 

We therefore selected our primary sample of dust- 
obscured galaxy candidates from the Bootes Field of the 
NOAO Deep Wide-Field survey based solely on the fol- 
lowing two criteria: 



> 0.3 mJy 
(R- [24]) > 14(Vega) mag 



(1) 
(2) 



There are 2603 sources satisfying these criteria in the 
Bootes field. Of these, 703 are fainter than the 50% com- 
pleteness limit of our optical imaging, and the (R— [24]) 
values are therefore lower limits. We refer to this set 
of 2603 galaxies as our "total sample", and use it for 
estimating the infrared luminosity density contributed 
by this population. This large number of sources corre- 
sponds to an average surface density of w 0.089 sources 
arcmin~ 2 . 

As mentioned above, a subset of the Bootes Field has 
been mapped by IRAC (Eisenhardt et al. 2004). The 
total area of the combined Bw, R, I, IRAC 3.6/xm and 
4.5/im, and MIPS 24/im observations is (after addition- 
ally masking out bright stars and bad regions) 6.76 deg 2 . 

14 We define the 24/tm Vega magnitude as [24] = 
-2.5 logio (-F24,im/7.29Jy). F 24Mm = lmjy (0.3mjy) corresponds 
to a Vega magnitude of 9.66 (10.96). 




Fig. 1. — Top panel: (R- [24]) color vs. 24/im magnitude distribu- 
tion for all sources in the NDWFS Bootes Field. The bottom and 
top abscissae show the 24/tm magnitude and flux density respec- 
tively, and the left and right ordinates show the color in magnitudes 
and the i^^m/ 1 Fr flux density ratio respectively. The objects with 
(-R-[24])>14 lie above the horizontal line. Bottom panel: The ex- 
pected variation of color with redshift of two local ULIRGs and the 
nearby starburst galaxy M82. (i?-[24])>14 sources have no local 
counterparts. 
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Fig. 2. — The differential sky density (top panel) and fraction 
(bottom panel) of 24/im sources with colors (R — [24]) > 14 mag 
as a function of 24/im flux density. The fraction rises steadily as a 
function of decreasing 24/im flux density beyond i 7 24pm ~ lmjy. 

All of the galaxies selected by criteria (1) and (2) above 
that lie within this region are detected in at least the 
3.6/im and 4.5/im IRAC bands. We refer to this result- 
ing sample of 1,882 sources as the "photometric sample" . 
This latter sample is used below in the discussion of the 
IRAC color-color diagrams and for constructing robust 
spectral energy distributions. 

2.3. Spectroscopic Data 

Members of this extremely red 24/im source popula- 
tion were targeted for spectroscopic observations using 
the Infrared Spectrometer on Spitzer (IRS; Houck et al. 
2004), the Low Resolution Imaging Spectrometer (LRIS; 
Oke et al. 1995), the Deep Imaging Multi-Object Spec- 
trograph (DEIMOS; Faber et al. 2002), and the Near 
Infrared Spectrometer (NIRSPEC; McLean et al. 1998) 
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on the telescopes of the W. M. Keck Observatory, and 
the Near-Infrared Imaging Spectrograph (NIRI) on the 
Gemini-N telescope. The Bootes field is also the tar- 
get of the AGN and Galaxy Evolution Survey (AGES; 
Kochanek et al. 2008 in preparation), which has mea- 
sured over 20,000 redshifts for optically bright galaxies 
(7<20) and AGN (I < 21.5) using Hectospec (Fabricant 
et al. 2005) on the MMT. The IRS observations focused 
primarily on the mid-IR bright sources which were un- 
detected at optical wavelengths, targeting objects with 
-p24/im~0.75 mJy; preliminary results have been pre- 
sented by Houck et al. (2005) and Weedman et al. (2006). 
Redshifts were measured from the mid-IR spectra by fit- 
ting either the strong (rest-frame) 9.7/zm silicate absorp- 
tion feature, or the 7.7/mi PAH emission features. Typ- 
ical uncertainties for the IRS redshifts are Az< ± 0.2. 
We refer the reader to the Houck et al. (2005); Weedman 
et al. (2006) and Desai et al. (2008, in prep) papers for 
details. 

The Keck optical observations targeted members of the 
(R — [24]) > 14 population as well as randomly selected 
24fim sources down to the 24^m flux density limit of 
0.3 mJy. The details of the Keck spectroscopic observa- 
tions will be presented elsewhere (Soifer et al. 2007, in 
preparation) and the resulting redshift distribution of the 
overall -F24^m > 0.3 mJy source population is presented 
in Desai et al. (2008, ApJ, submitted). Redshifts of 
the sample discussed in this paper were derived from the 
optical spectroscopy by (typically) identifying either the 
[OII]AA3726,3729 doublet in the far red (for the z<1.6 
galaxies) or the Lya line in the blue for the z> 1.8 ob- 
jects. The near-infrared spectroscopic redshifts are based 
on the detections of the Ha+ [Nil] or [OIII] emission lines 
(see Brand et al. 2007, for examples). Typical uncertain- 
ties for the optical and near-IR measured redshifts are 
Az< 0.002. 

Based on these various observations (but primarily on 
the Spitzer and Keck efforts), we have now measured 
a total of 86 spectroscopic redshifts for the sample of 
candidate dust-obscured galaxies (i.e., i<24p m > 0.3 mJy 
and (R - [24]) > 14) in the Bootes field. Roughly half 
the redshifts (41/86) are the result of Spitzer/IRS obser- 
vations, and half come from ground-based optical (38) 
or near-IR (7) spectroscopy. There are only 5 galaxies 
which have redshift measurements from both Spitzer/IRS 
and Keck. These 5 have a mean redshift difference of 
Az = ZKcck — ~ —0.04 (with an uncertainty on the 
mean of ±0.15) and a standard deviation of <ta 2 ~ 0.3. 
Hence, we do not detect any statistically significant sys- 
tematic redshift offset for the IRS redshifts, and adopt 
them as the true redshifts of the galaxies for the purposes 
of this paper. 

We also compared the spectroscopically measured red- 
shifts with photometric redshift estimates from various 
template fitting algorithms, and found a very large scat- 
ter between photometric redshift estimates and the spec- 
troscopic measurements; photometric redshifts are not 
reliable for this extreme population. Given the uncer- 
tainties inherent in fitting templates to spectral energy 
distributions dominated by power laws and only weak 
features, we do not attempt to estimate photometric 
redshifts (cf. Perez- Gonzalez et al. 2005; Polletta et al. 
2006), and instead use the spectroscopically measured 
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Fig. 3. — Same as figure 1, but restricted to sources with colors 
(R— [24]) > 14. The dots represent the entire sample and the solid 
(blue) circles represent the galaxies with measured spectroscopic 
redshifts. The upward arrows show sources with upper limits on 
their ij-band magnitudes. The spectroscopic redshifts sample the 
entire range of color and magnitude. 

redshift distribution directly. 

The galaxies for which we have measured redshifts 
sample the full range of F 2 4 l _ lm and (R — [24]) color (Fig- 
ure 3). For the objects with measured (R — [24]) colors 
(i.e., excluding the galaxies with only upper limits), the 
color distribution of galaxies with redshifts is a reason- 
able sampling of the full color distribution of the sample 
(i.e., a two-sided Kolmogorov-Smirnov test finds that the 
color distributions are similar at the 38% level). The 
sampling of the 24/xm flux density distribution is much 
poorer. This is because a significant fraction of the source 
redshifts were measured with Spitzer/IRS, and thus are 
restricted to the galaxies with F 2 4^ m ^0.5 mJy. Nev- 
ertheless, the spectroscopic sampling does span the full 
range in both color and 24/im flux density, and we as- 
sume for the purposes of this paper that the resulting 
redshift distribution is a fair representation of the over- 
all population. 

3. RESULTS 

In this section we discuss the colors and overall spectral 
energy distributions (hereafter SEDs) of the R— [24] > 14 
galaxy population. We also present the redshift distribu- 
tion, and under the assumption that the measured distri- 
bution applies to the entire population, we estimate the 
rest- frame UV, mid-IR and IR luminosity distributions. 
Finally, we discuss the X-ray emission from the sample, 
with the aim of understanding whether AGN play a role 
in the energetics. 

3.1. Colors and SEDs 
3.1.1. Mid-Infrared Colors 

Several recent studies have demonstrated that the 
mid-infrared colors of galaxies can provide an effective 
discriminant between 'normal' galaxies and unobscured 
AGN (e.g., Lacy et al. 2004; Stern et al. 2005). The top 
panel of Figure 4 shows the distribution of mid-infrared 
color for all the 24/im sources with i 7 24 Aim > 0.3 mJy 
in the Bootes NDWFS field. The color distribution is 
similar to that seen by Stern et al. (2005), showing a 
sequence of low-redshift star-forming galaxies at 'blue' 
[3.6] — [4.5] colors, and a second sequence at redder col- 
ors that largely overlaps the 'AGN wedge' (shown by the 
dashed line; see Stern et al. 2005). 
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Fig. 4. — (a) The mid-infrared color distribution of all 24/im 
sources in the NDWFS Bootes field with flux densities ^24^01 > 
0.3 mjy, shown as a grey scale. This greyscale is repeated in all 
panels. The dashed lines demarcate the region occupied by bright 
AGN as defined by Stern et al. (2005). (b) The mid-infrared color 
distribution of the DOGs ((R — [24]) > 14 sources; shown as solid 
circles) with 24/xm flux densities F^A^m > 1 mJy- (c) Same as (b), 
but for DOGs with 0.6 < i*24um < 1 mjy. (d) Same as (b), but for 
DOGs with 0.3 < F24/jm < 0.6 mjy. The DOGs span a large range 
of [5.8]-[8.0] color, and the color is a function of the observed 24/im 
flux density, with the fainter sources having bluer colors. Notice 
that there is a continuum of sources across the "AGN boundary" 
of Stern et al. (2005) at ([5.8]-[8.0])=0.6. 

The mid-infrared colors of the sources with (R— [24]) > 
14 are shown in the lower panels of figure 4. The 
(R - [24]) > 14 population lies at redder [3.6]- [4.5] col- 
ors (with the majority of objects having [3.6] — [4.5]>0.5), 
suggesting that this population lies at high rcdshift 
(z>l) or are AGN (Stern et al. 2005). Their positions 
in the color-color diagram also overlap the 'AGN wedge', 
but spill across the boundary of the wedge toward bluer 
[5.8] — [8.0] colors. The figure also demonstrates that the 
location of the population in the color-color diagram is a 
function of -pM^m, with the apparently brighter sources 
(with F 2 4 flm ^ lmJy) lying toward redder [5.8] — [8.0] col- 
ors, and the fainter sources (with i 7 24^m J; 0.6mJy) hav- 
ing bluer [5.8] — [8.0] colors. The sequence appears to be 
continuous across the boundary of the 'AGN wedge' at 
([5.8] — [8.0])=0.6, with no obvious break in the popula- 
tion. The median values of the ([3.6]-[4.5],[5.8]-[8.0]) 
colors in the three successive bins of increasing 24/im 
flux density are (0.70,0.64), (0.79,1.07) and (0.95,1.40) 
mag. A two-sided Kolmogorov-Smirnov (hereafter K-S) 
test rules out the hypothesis that the color distribution 
in one of the three flux density bins is drawn from any 
other one at a very high level of significance. 

3.1.2. Overall Spectral Energy Distributions 

The SEDs of galaxies can provide important insights 
into their nature, i.e., the components that dominate the 
energetics at different wavelengths. The overall SEDs of 
the sources with (R — [24]) > 14 span a wide range. At 
one extreme, the SEDs are fairly monotonic throughout 
the mid-IR bands, and are fairly well represented by a 
power-law, suggesting a broad range of (warm) dust tem- 
peratures similar to that observed in AGN-dominated 
sources. At the other extreme, the SEDs show a dis- 
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Fig. 5. — Examples of the mid-infrared spectral energy distri- 
butions observed for the [R — [24]) > 14 sample. The SEDs range 
from being power-law dominated (near the top) to exhibiting a 
distinct "bump" in their SED between 3— 10/tm (near the bottom). 
The "bump" is attributed to starlight in the galaxies, peaking at 
a rest-frame wavelength of ssl.6/tm. 

tinct "bump" , which we interpret as being due to the 
presence of a stellar continuum peaking at a rest-frame 
wavelength of 1.6/im. Examples of the two classes of 
SEDs arc shown in figure 5. There is a continuum of 
SEDs between these two extremes, and it is difficult to 
define an objective criterion that cleanly discriminates 
between these two classes, perhaps suggesting that the 
SEDs in a large fraction of the population are character- 
ized by both AGN and star-formation activity. 

In order to establish a more rigorous definition (more 
for statistical than physical purposes), we define the fol- 
lowing objective means of classifying a source SED as 
being power-law dominated. We performed two power- 
law fits to the mid-infrared flux measurements of every 
source, the first to just the IRAC measurements (i.e., the 
3.6/im to 8.0/im data) and the second to the combined 
IRAC and MIPS 24^im data. We then examined the 
power-law indices (oirac and aiRAC+24; Fu ^ and 
X 2 values of the resulting fits. An SED was classified as a 
'power-law' source if the power-law index Chirac > 3, or 
if ckiRAc > 1 and Xirac ^ 1-5 , or if oirac > o:t rac+24 

and |aiRAc - ai RA c+24| > 2^ ^irac + <J mAC+24' wh ere 
f irac j cirac+24 are the formal fit uncertainties on the 
a's. Although this selection appears complex, it does 
separate the sources which show monotonic SEDs from 
the ones which show "bumps" or are otherwise different 
(due to being too noisy, for example). 

Using this SED discriminant, we find that the spec- 
tral energy distributions are loosely correlated with 24/im 
flux density, with the bright -F^m sources being predom- 
inantly power-law sources (presumably with AGN emis- 
sion dominating the mid-infrared spectral region) and 
the fainter ^24^111 sources exhibiting clearer 'bumps' in 
their SEDs (Figure 6). 

Our SED discriminant is constructed solely to capture 
this qualitative trend in the observed SEDs, and not to be 
rigorously used to define unambiguous samples of 'bump' 
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Fig. 6. — The fraction of power-law sources as a function of 24/im 
flux density for the (R — [24]) > 14 sample of galaxies. See text for 
details. 

or 'power-law' sources. Different thresholds in \ 2 or a 
in the criteria above change the fractions of power-law 
sources, but still reproduce the basic trend. 

The range of SEDs in our sample is similar to that 
observed by Polletta et al. (2006). Polletta et al. (2006) 
define several classes based on the shape of the SED, and 
identify "class I" and "class II" sources as their primary 
sample of obscured infrared AGN. Our (R — [24]) > 14 
criterion would select all of the "class I" sources and 36% 
of the "class II" sources in the Polletta et al. (2006) study. 

3.2. Redshifts 

The redshift histogram of all 86 galaxies is shown in 
Figure 7. The distribution is broad, and can be modeled 
as a Gaussian centered at z = 1.99 ± 0.05 with a disper- 
sion a z = 0.45 ±0.05. The overall redshift distribution is 
very similar to that of galaxy populations selected at sub- 
millimeter wavelengths (the Sub-mm Galaxies, or SMGs; 
Chapman et al. 2005). 

It is noteworthy that the redshift distribution of galax- 
ies with 14 < (R— [24]) < 15 differs statistically from that 
of redder galaxies with (R— [24]) > 15. Both a Student's 
T-test and a Wilcoxon Rank-Sum test rule out the hy- 
pothesis that the two subsamples have similar means at 
>99.9% significance, and a two-sided K-S test finds that 
the two distributions differ at the «99.6% level of signif- 
icance. The bulk of the redshifts for the reddest galaxies 
(i.e., R — [24] > 15) were obtained with Spitzer/IRS, 
but since the distributions of the IRS-determined red- 
shifts also show a similar dependence on R — [24] color, 
this observational bias does not appear to be the source 
of the difference (see the red shaded histograms in Fig- 
ure 7). Unfortunately the redshift distribution is too 
broad, and the spectroscopic redshifts too few, to trust 
any correlation between the (R — [24]) color and the 
redshift derived from the present data, but it does ap- 
pear that our color selection criterion effectively selects 
against lower redshift sources. Indeed, applying the color 
criterion (R— [24]) > 14(15) to the sample of Yan et al. 
(2007) effectively selects only the sources with redshifts 
z > 1(1.6). 

There is no obvious dependence of redshift on mid- 
infrared color or 24/im flux density within the DOG pop- 
ulation. In particular, we find no significant correlation 
between the [3.6] -[4.5] or [5.8] -[8.0] color with redshift 
within our sample (Figure 8). A two-sided K-S test finds 




Redshift 

Fig. 7. — Top: The observed redshift distribution of all 86 
(i? — [24]) > 14 galaxies. The dashed line in the top panel shows 
the best fit Gaussian approximation to the observed redshift distri- 
bution, with z = 1.99 ± 0.05, a = 0.45 ± 0.05. Lower three panels: 
Redshift distributions in more restricted ranges of (R — [24]) color, 
as designated in the panel legends. The shaded histogram in each 
of the three lower panels shows the subset of redshifts measured 
from Spitzer/IRS observations. 

no significant difference in the redshift distributions for 
galaxies brighter than or fainter than -FiM/xm = 1 rnJy. 
There is some evidence that the small subset of 22 
galaxies with F24 flm < 0.6mJy does show a narrower 
redshift distribution centered at slightly lower redshifts 
(z(F24 pm < 0.6mJy) sa 1.86). This trend is also rein- 
forced when considering the redshift distribution of all 
24/im sources (irrespective of color). Desai et al. (2007) 
present the redshift distribution of a sample of galax- 
ies selected with 24/im flux densities i*24/un > 0.3mJy, 
and demonstrate that there is general trend of increasing 
mean redshift with (R— [24]) color. However, our current 
sample of 86 redshifts for (R— [24] ) > 14 sources is a small 
subset, and without a larger spectroscopic sample, we 
are loathe to adopt anything other than the simplest as- 
sumptions regarding the redshift distribution. Since the 
current set of spectroscopic redshifts, although obtained 
in diverse ways, span the full range of optical magnitude 
and 24/im flux density of the overall DOG sample (see 
Figure 3), we adopt the Gaussian fit to the observed red- 
shift distribution as the overall redshift distribution of 
the sample. 

The redshift distribution implies that the observed 
R — [24] selection corresponds (approximately) to a rest- 
frame 8/xm to 2200A flux density ratio. In order to search 
for local counterparts of these galaxies, we examined the 
Spitzer/WAC to GALEX/NVV flux density ratios of all 
550 low-rcdshift (z < 0.1) galaxies in Bootes with spec- 
troscopic redshifts from the AGES survey (Kochanek et 
al., in prep); none were found to have flux density ratios 
similar to those in our DOG sample (S. Salim, personal 
communication) . 

3.3. Rest-Frame Luminosities 
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Fig. 8. — The mid-infrared [3.6/im] — [4.5/im] (top panel) and 
[5.8/nm] — [8.0/^m] (bottom panel) colors as a function of redshift 
for DOGs with measured spectroscopic rcdshifts and > 3cr mea- 
surements of flux densities in the IRAC bands. The dashed line in 
the lower panel represents the lower color cut used by Stern et al. 
(2005) to select AGN. There is no strong correlation between color 
and redshift for the R — [24] > 14 sample. We interpret the large 
spread in color at a given redshift as representing a range in the 
relative contributions of reprocessed AGN and stellar luminosity to 
the observed mid-infrared bands, with Lmir(AG N) / LMm{star s) 
increasing toward redder (i.e., warmer) colors (cf. Figure 11 of Maz- 
zarella et al. 1992). Rather than a clean distinction between AGN 
and star-forming galaxies, there is a continuum of sources across 
the range of color. 

3.3.1. 8[im Luminosities 

The rest- frame 8/im luminosity, Lg, of the DOGs is 
estimated using 



4< 



Lg = vL^g^ — ^8jjm j-— — Fs^m 



(3) 



+ z 

^8ftm ~ F 2 4nm ^ ^ ^ , (4) 

where a = 2.096 log(F 2 4 Aim /F 8Aim ) is the spectral index 
between 24/xm and 8^m, dj, is the luminosity distance, 
and vg^m is the frequency corresponding to 8/xm. 

The 8^m luminosities for the DOGs with measured 
spectroscopic redshifts are shown in figure 9. Our esti- 
mates assume that the spectral shape between (observed 
frame) 8^m and 24/im is well matched by a power-law; 
this may be true if the mid-infrared emission from the 
galaxies is dominated by a featureless AGN, but not if 
the spectrum is strongly affected by PAH emission or 
silicate absorption. 

The derived Lg values are large, with most of the 
galaxies being more luminous than lO n L and half with 
Lg > 10 12 L Q ; by comparison, the 8/xm luminosity func- 
tion derived by Reddy et al. (2007) for z 2 UV-selected 
star-forming galaxies only extends to ~ 10 n LQ. The 
median 24/xm flux_density of the sources with spectro- 
scopic redshifts is F 2 4^ m — 1.04mJy, in comparison with 
the entire sample, which has F 2 4^ m ~ 0.4. Hence, if the 
overall population of galaxies with R — [24] > 14 has the 
same redshift distribution as the subset of galaxies with 
spectroscopic redshifts, then the median rest-frame 8/im 
luminosity of the overall population is Lg 4 x 10 n L Q . 

3.3.2. Far-Infrared Luminosities 
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Fig. 9. — Distribution of rest-frame 8/im luminosity (fL ly (8fj,m)) 
of all 86 sources in the (R — [24]) > 14 samples with spectroscopi- 
cally measured redshifts. The top abscissa shows the corresponding 
Lir, estimated using the conversion of Caputi ct al. (2007). The 
bulk of the sources are clearly ULIRGs. (Our adopted conversion 
between Lg and -Ltr is only accurate to within a factor of three; 
to use a more [less] conservative conversion of Ltr, = 5[15]Ls, add 
~ 0.70 [1.18] to the lower abscissa.) 

The distribution of far-infrared luminosities (Iir) for 
the DOGs with spectroscopically measured redshifts is 
shown in Figure 9 (see upper abscissa). Nearly all (80 
/ 86, or 93%) of the sources with spectroscopic redshifts 
are ULIRGs, with Ltr > 10 12 L Q , and more than half of 



these have Lir > 10 13 L Q . 

Converting from observed 24^m flux densities to rest- 
frame far-infrared or bolomctric luminosities is notori- 
ously uncertain and has a checkered history (e.g., see 
§6 of Caputi et al. 2007, for a recent summary). The 
conversion factors vary significantly depending upon the 
template adopted for the extrapolation, ranging from 
L m «5L 8 for Mrk 231 to L m w 27.6 Lg for Arp 220. 
(Here, we adopt the Sanders & Mirabel (1996) defini- 
tion of Lir = L(8/im — 1mm).) As pointed out by vari- 
ous authors (e.g., Reddy et al. 2007; Caputi et al. 2007), 
many studies of the star-forming galaxy population have 
assumed the higher values, especially for the most 8/xm- 
luminous galaxies. In constructing the upper abscissa 
of Figure 9 we used the relation determined by Caputi 
et al. (2007) for the conversion, i.e., L m = 1.91 L\ m , 
where Lg = vL^g^ in units of Lq. However, this re- 
lation was derived by fitting [L 8 7^ir] measurements of 
93 low-redshift (0 < z<0.6) galaxies (from the Spitzer 
First Look Survey), very few of which have Lg > 10 n Lo 
or Lir > 10 12 . Hence the extrapolation of this relation 
to higher redshift and higher luminosity galaxies is ques- 
tionable. 

Therefore, we also quote a range of Ltr, values deter- 
mined from Ltr = 10^ 5 Lg. This range encompasses the 
conversion estimate based on Mrk 231 at the low end, 
and that adopted by Reddy et al. (2006) for z <~ 2 galax- 
ies at the high end. Since the mid-infrared emission from 
the galaxies in our sample appears to be dominated by 
AGN at the bright end and star-formation at the faint 
end, the large range we adopt should encompass the true 
values for Lir. Clearly, measurements of the rest-frame 
far-infrared spectral energy distributions of these galax- 
ies are necessary before more precise estimates of Lir 
can be made (cf. Tyler et al. 2008, ApJ submitted). The 
high far-infrared luminosities for DOGs hold even if we 
adopt the more conservative scaling based on Mrk 231: 
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Fig. 10. — Distribution of i?-band (Vega) magnitudes for all 
the (R — [24]) > 14 galaxies (shown as the solid-line histogram) 
in the Bootes Field, compared with the distribution for the z 2 
"BX" UV-selected star-forming galaxies from Reddy et al. (2007) 
(shown as the dashed-line histogram). The subset of our sample 
with only upper limits in the i?-band are represented by the shaded 
histogram. The vertical dotted line represents the magnitude limit 
of the "BX" sample. The (R — [24]) > 14 galaxy sample lies to 
fainter magnitudes than the "BX" sample, with half the sample 
fainter than the magnitude limit of the "BX" sample. The his- 
tograms shown are not normalized for the survey areas. 

if Lib, = 5 L$, then 75 / 86, or 87% of the sources with 
spectroscopic redshifts are ULIRGs. 

3.3.3. Rest-Frame UV Luminosities 

By selection, the galaxies in our sample are faint at 
optical wavelengths. Figure 10 shows the distribution of 
optical apparent magnitudes for the (R— [24]) > 14 sam- 
ple. Most of the galaxies are fainter than i?=24.5 mag, 
and «27% of the galaxies have only upper limits for the 
i?-band magnitudes. For comparison, we also plot the 
distribution of i?-band magnitudes for the z ps 2 UV- 
selected star-forming galaxies from Reddy et al. (2007) 
(selected using their "BX" selection criteria). The ma- 
jority of galaxies in our sample lie at fainter i?-band mag- 
nitudes than the bulk of the "BX" galaxy sample, with 
roughly half having magnitudes fainter than the magni- 
tude limit of the "BX" selection. 

Figure 11 shows the distribution of absolute magni- 
tudes (at 2200A) for the subset of our sample with 
measured spectroscopic redshifts. Here, the abso- 
lute magnitude is estimated using M ° = Rab — 

° ° 2200A 

51ogio(di/10parsec) — 2.51ogio(l + z). The median UV 
luminosity of our population is M ° ps —19.7. The 

J 1 ^ 2200A 

distribution of UV luminosities for the "BX" sample of 
Reddy et al. (2007) is also shown for comparison. Nearly 
all the (R — [24]) > 14 galaxies are fainter than L* 
(M* ps -21.0AB mag, as determined from the "BX" 
sample of Reddy et al. 2007) , and many are approaching 
the UV rest- frame luminosities of local dwarf galaxies. 
Traditionally, low UV luminosity galaxies are considered 
to be less bolometrically luminous as well, and thought 
to contribute less to the overall IR luminosity budget. In 
contrast, the galaxies in our sample have low rest-frame 
UV luminosities but are nevertheless very bolometrically 
luminous and belong at the top end of the galaxy lu- 
minosity function. This underscores the importance of 
a multiwavelength approach in selecting and studying 
galaxy populations at high redshift. 
We have ignored the fc-corrections in our conversion 



Fig. 11. — The solid-line histogram shows the distribution 
of 2200A rest-frame absolute magnitudes (in AB units) of the 
(R— [24]) > 14 galaxies with spectroscopically measured redshifts. 
The subset of our sample with only upper limits in the i?-band are 
represented by the shaded histogram. The absolute magnitudes 
of the z £3 2 "BX" UV-selected star-forming galaxies is shown 
by the dashed line histogram, and the vertical dotted line repre- 
sents the magnitude of an L* galaxy. Nearly all the galaxies in 
the (R — [24]) > 14 sample lie below L* , and reach very faint 
rest-frame UV absolute magnitudes despite their remarkable mid- 
infrared luminosities. The histograms shown are not normalized 
for the survey areas. 

to absolute magnitude above, and assumed that all our 
sample galaxies lie at z — 2. This is defensible for 
the flat-spectrum galaxies in the "BX" sample of Reddy 
et al. (2007), but perhaps less so for our galaxy sample, 
whose members have very red spectral energy distribu- 
tions. Nevertheless, we have adopted this conversion for 
simplicity, and this choice makes little difference to the 
overall comparison of the two samples. 

3.4. X-Ray Properties 

Observations at X-ray wavelengths are well suited to 
identifying AGN over a range of redshift. The Bootes 
Field of the NDWFS has been mapped in its entirety by 
the Chandra X-ray Observatory (i.e., the XBootes Sur- 
vey; Kenter et al. 2005; Murray et al. 2005; Brand et al. 
2006). In order to determine whether the population of 
galaxies with extreme R — [24] colors is dominated by 
AGN or star-forming galaxies, we examined their X-ray 
emission. 

At a redshift z = 2, the Chandra soft (0.5—2 keV) 
and hard (2—7 keV) bands correspond to rest-frame en- 
ergies of 1.5-6 keV and 6-21 keV respectively. 2,101 
of the galaxies in our sample lie within 10' of the cen- 
ter of a Chandra X-ray Observatory ACIS observation 
in the XBootes Survey. Of these galaxies, 59 have 
X-ray detections with >4 counts in the XBootes sur- 
vey; this corresponds to an X-ray flux of -Fb.5-7keV ~ 
7.8 x 10 -15 erg cm -2 s _1 and, at z — 2, a luminosity 
of L x = 2 x 10 44 erg s" 1 (Brand et al. 2006). This lu- 
minosity is more typical of AGN than of star-forming 
systems (e.g., Colbert et al. 2004). As shown in the up- 
per panel of figure 12, the X-ray detection rate of the 
galaxies in our sample increases with 24/zm flux den- 
sity, from ps 1% at F 2 4 Mm ps 0.3 mJy to > 10% at 

^24/im <^ 1 mJy. This detection rate is similar to that of 
the overall 24/im source population (shown in the lower 
panel of Figure 12) at bright F^m, but a factor of two 
lower at the faint end, suggesting that the AGN frac- 
tion within the (R — [24]) > 14 sample is an increasing 
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Fig. 12. — Top panel: The fraction of galaxies in the (R — [24]) > 
14 sample detected as X-ray sources in the Chandra XBootes survey 
(Murray et al. 2005; Kenter et al. 2005; Brand et al. 2006), shown as 
a function of 24/im flux density. Sources which are brighter at 24/im 
have a higher X-ray identification rate. Bottom panel: Fraction of 
all 24/im sources detected as X-ray sources in the XBootes Survey. 
The X-ray identification fraction in the (R — [24]) > 14 sample is 
similar to that of all 24/im sources at the bright end, but smaller 
at faint flux densities F24 Mm ~ 0.3 — 0.6 mJy. 

function of 24^m flux density. 

The X-ray data suggest that even the fainter galax- 
ies may contain a significant AGN component, despite 
displaying oft-used discriminants of star-forming galax- 
ies (e.g., stellar "bumps" in their SEDs, colors typically 
attributed to star-forming galaxies, low values of the 
24/im to 8/im flux ratio, and mid-infrared spectra ex- 
hibiting PAH emission). Although there may be AGN 
in the fainter sources, it is unclear from the current data 
whether the AGN plays a significant role in the mid/far- 
infrared luminosity of the fainter i 7 24/xm galaxies, but if 
it does, then the global star formation rate density de- 
rived by including the full contributions of their mid/far- 
infrared luminosity may be overestimates. 

4. DISCUSSION 
4.1. The Nature of the Extreme Red Population 

What arc the extreme red objects that become an 
increasing fraction of the population at fainter 2Afim 
flux densities? The spectroscopically measured red- 
shifts for this population suggest prodigious luminositcs 
(10 12 ~ 14 L Q ) anc | j-gddgj. colors than any local ULIRGs. 
We therefore suggest that the rise in this red fraction 
is predominantly an evolutionary effect due to the ap- 
pearance of a more enshrouded population of luminous 
objects at higher rcdshift. The space density of this pop- 
ulation is low; if the comoving space density remains con- 
stant with time, they may evolve into the fairly luminous 
local ellipticals. We speculate that we have uncovered a 
population of young galaxies undergoing very rapid bulge 
formation and / or AGN accretion, that evolve into fairly 
luminous galaxies (<~ AL*) in the present-day universe. 

The optical-mid-IR SEDs of this population show a 
range of color. However, by selection, these objects 
are extreme in their optical-to-mid-IR properties and, 
as demonstrated by Figure lb, have colors redder than 



any local dusty galaxy or ULIRG. Since we have only 
limited spectral information, the dust reddening can- 
not be determined unambiguously. However, the exist- 
ing mid-infrared spectroscopy shows that these galaxies 
have spectra that rise into the mid-infrared, suggesting 
emission from warm dust, and exhibit cither deep sili- 
cate absorption (e.g., Houck et al. 2005; Weedman et al. 
2006) or PAH emission (Yan et al. 2007, ; Desai et al. 
2008, in prep). Near- infrared spectroscopy by our group 
for a handful of galaxies reveals, where measurable, large 
Balmer Ha/H/3 decrements in both the narrow and broad 
components of the Balmer lines (Brand et al. 2007). This 
suggests that the galaxies contain significant quantities 
of dust distributed on scales of, at the very least, the 
extended (i.e., kiloparsec scale) narrow-line emitting re- 
gion. Since these galaxies lie at high redshift and are 
therefore very luminous, they must be undergoing sig- 
nificant AGN accretion and / or star-formation activ- 
ity. This would suggest that the large rest-frame near- 
infrared to ultraviolet luminosity ratios result from heavy 
extinction of the most luminous emitting regions. 

It is difficult to estimate the extinction in these galaxies 
from the limited photometric information, but we can de- 
rive some general constraints by considering the extinc- 
tion necessary to create such red SEDs given a young stel- 
lar population. Figure 13 shows the U-band extinction 
(in magnitudes) necessary to create an (R — [24]) = 14 
color for template galaxies with ages less than 1 Gyr. 
This is based on the assumption of an extinction law de- 
rived using the functional form of Calzetti et al. (2000) 
for A < 2.2^m and values of the extinction as tabulated 
in Table 4 of Draine (2003) at longer wavelengths. Al- 
though the assumption of a monolithic dust screen is 
unrealistic, it does provide a lower limit on the extinc- 
tion required in these galaxies. In addition, the fact that 
the entire galaxy is so underluminous at rest-frame UV 
wavelengths suggests that either the stellar population 
in these systems is old (i.e., <J 1 Gyr) and intrinsically 
faint at UV wavelengths, or that the dust is distributed 
on large enough scales to shroud the stellar component 
as well. If these galaxies are young and dominated by a 
star-forming phase, then they are enveloped on kilopar- 
sec scales by prodigious amounts of dust. 

Given the redshifts (z <~ 2) and spectral energy distri- 
butions of these galaxies, it is very unlikely that they are 
dominated by old stellar populations, and we therefore 
hypothesize that these objects are dust-obscured galax- 
ies (DOGs), where most of the luminous regions of the 
galaxy are shrouded in dust. Such a scenario would sug- 
gest an overall SED for these galaxies containing three 
(perhaps related) components: a reddened stellar com- 
ponent peaking at a rest-frame wavelength of 1.6^m, cold 
dust enshrouding the star-forming regions rcradiating 
the absorbed UV emission at long wavelengths (corre- 
sponding to, say, a rest-frame wavelength of 60^m), and 
warm dust enshrouding the AGN dominating the emis- 
sion in the rest-frame mid- infrared (i.e., at rest-frame 
wavelengths of 5-20/im). 

If the luminosity from the sources selected by (R — 
[24]) > 14 is largely the result of star formation activ- 
ity, then the implied star-formation rates are very large 
1700M Q /yr for an L m = 10 13 L© object; Kennicutt 
1998). Such a violent starburst is unlikely to be confined 
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Fig. 13. — The minimum extinction required to redden the spec- 
tral energy distributions of young stellar populations and AGN to 
produce an optical-to- mid-infrared color of (R — [24]) = 14. The 
template SEDs are models from the Bruzual & Chariot (2003) pop- 
ulation synthesis library with ages of 5 Myr (solid line), 25 Myr 
(dotted line), 100 Myr (short-dashed line) and 900 Myr (dot-dashed 
line), and the median QSO template (long-dashed red line) from 
Elvis et al. (1994). The reddening law used is a combination of that 
from Calzetti et al. (2000) and longer wavelength estimates from 
Draine (2003), and assumes the unrealistic case of a dust screen in 
front of the emitting source in order to derive a firm lower limit on 
Ay. 

to a small region, and the (observed) very large infrared- 
to-ultraviolet luminosity ratio must result from dust dis- 
tributed on very large (kiloparsec) scales. In addition, 
star-formation rates this large cannot be sustained for 
long periods, suggesting that these galaxies are viewed 
during a short-lived phase during which the bulk of the 
stars is produced. 

If, instead, the sources are dominated by AGN, then 
the dust could be mostly confined to the nuclear regions 
of the galaxy. However, based on the measurement of 
large Balmer decrements in both broad- and narrow-line 
components in a few galaxies, Brand et al. (2007) ar- 
gue that the dust distribution must be patchy and over 
large volumes. Even if the galaxy luminosity is domi- 
nated by AGN activity, this phase of evolution must be 
fairly short lived. The bolometric luminosity of an AGN 
may be written (in terms of its Eddington Luminosity) as 
L^ N ee eL Edd » 0.33 x lO 13 (e/O.l)(Af BH /lO 9 Mp)L , 
or in terms of its accretion rate as L&£[ = rjMc 2 rs 
0.15 x 10 13 (?7/0.1)(Af/lM Q yr^ 1 )L Q . Large AGN lumi- 
nosities therefore generally imply high accretion rates 
onto massive black holes. The present data on the DOGs 
do not allow any useful constraints on either the radiative 
efficiency, the accretion efficiency, or the black hole mass, 
but their luminosities of 10 13 Lq suggest that they must 
be radiating at close to the Eddington luminosity even if 
they harbor a very massive black hole. Even at high ac- 
cretion efficiencies (i.e., r\ = 0.1), the timescales to build 
up a 10 9 Af© black hole would be short («0.15 Gyr). A 
flux density limited sample of galaxies (such as ours) will 
preferentially identify the galaxies in their most luminous 
phase, suggesting that we may be witnessing the last 
throes of the AGN formation in these obscured galaxies. 

What is the relevance of this rather extreme popula- 
tion to the overall galaxy population at z « 2? In the fol- 
lowing subsections, we estimate the space density of the 
DOGs and their contribution to the infrared luminosity 
density in an attempt to address this question. In con- 



trast to previous work, we do not discriminate between 
AGN and star-forming galaxies within our sample, but 
estimate their joint contribution to the space densities 
and contribution to the infrared luminosity density. 

4.2. Space Density of DOGs 

Based on the observed redshift distribution, we now 
estimate the space density of this population. We assume 
that the comoving volume sampled by the (R— [24]) > 14 
and ^24^111 > 0.3 selection criteria is simply the fraction 
of the volume between 0.5 < z < 3.5 weighted by the 
redshift distribution. In other words, the effective co- 
moving volume sampled is given by: 



U c cff (DOGs) = ASlfg( 



z)dV c 



where Af2 is the solid angle covered by the survey, 



9(z) 



1 



=exp 



(z-z) 



2(7? 



(5) 



(6) 



<J Z V27T 

(with z — 1.98 and a z = 0.53) is the (normalized) Gaus- 
sian fit to the redshift distribution, and dV c is the comov- 
ing volume element per unit solid angle per unit red- 
shift interval (e.g., Peebles 1993; Hogg 1999). For the 
Bootes field, the volume sampled by the redshift distri- 
bution is roughly a third of the comoving volume between 
0.5 < z < 3.5. The resulting average space density for 
this population in the redshift range 0.5 < z < 3.5 is 
S D OG(^24 M m > 0.3 mJy) = (2.82±0.05) x 10" 5 /i? Mpc" 3 . 

(7) 

This space density of DOGs is comparable to that of 
highly UV-luminous z ~ 2 — 3 star-forming galaxies se- 
lected using the Lyman-break or near-infrared selection 
techniques (e.g., Reddy et al. 2007; van Dokkum et al. 
2006, and references therein). For example, the dust- 
obscured galaxy population has a space density compa- 
rable to that of Af A B(1700A) - 51og/i .7 ~ -22.2 UV- 
bright galaxies at z ~ 1.9 — 3.4, based on the lumi- 
nosity functions derived by Reddy et al. (2007). Al- 
though the space densities are comparable, it is impor- 
tant to note that each DOG contributes significantly 
more IR (and bolometric) luminosity per galaxy to 
the universe at z ~ 2. The population of "massive 
galaxies" selected as "Distant Red Galaxies" (DRGs) at 
near- infrared wavelengths has a space density of (2.2 ± 
0.6) x 10" 4 /i^ galaxies Mpc -3 (van Dokkum et al. 2006), 
roughly ten times more numerous than DOGs. As men- 
tioned above, the DRG and DOG populations also over- 
lap in their rest-frame UV luminosity distributions, and 
it is therefore possible that a fraction of the star-forming 
DRGs are selected by our (R — [24]) > 14 criterion. 

A recent measurement of the luminosity function by 
Wake et al. (2006) finds that luminous (i.e., > AL*) 
red galaxies in the local Universe have a space density 
of 2.5 x 10~ 5 /i 3 Mpc~ 3 , very comparable to that of 
the z w 2 DOGs. It is tempting to suggest that DOGs 
may therefore represent a common phase of evolution for 
most massive galaxies. Indeed, if this phase is short lived 
(i.e., if the duty cycle of the enshrouded, luminous phase 
when galaxies may be selected as DOGs is short), then 
the DOGs may be the progenitors of a large fraction of 
present-day > L* galaxies. 

More detailed comparisons will require better deter- 
minations of the redshift distribution as a function of 
apparent magnitude and 24/im flux density. 
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4.3. Infrared Luminosity Density of DOGs 

In order to estimate the infrared luminosity density 
contributed by the DOGs, we begin by assuming, as 
before, that all the DOGs with no measured spectro- 
scopic rcdshifts have a redshift probability distrbution 
described by the Gaussian approximation to the observed 
redshift distribution. For the DOGs with spectroscopic 
rcdshifts, we computed their L s values and converted 
them to Lir estimates based on the procedure described 
above. For the DOGs without measured redshifts, we 
assigned random rcdshifts (drawn from the probability 
distribution) to each object and computed their lumi- 
nosities. As in section §3.3.2, we computed luminosities 
using the conversion of Caputi et al. (2007), and derive 
a range based on the two extremes Lib, = [5 — 15] x L%. 
Since DOGs with fainter are also typically fainter 

in the IRAC bands, we adopt the mean a = 2.296 as the 
spectral index for all sources. 

The infrared luminosity density (IRLD) of the entire 
sample was estimated using 1000 Monte Carlo simula- 
tions of the redshift distribution of the galaxies. We find 
the total infrared luminosity density contributed by the 
DOG population to be 

log[IRLD DO G(F24 M m > 0.3mJy)] = 8.2281°;™, ( 8 ) 

where IRLD is in units of L© Mpc~ 3 and the limits are 
not ler uncertainties, but the likely range of the IRLD 
based on the uncertainties in the Lg, to Lir conver- 
sion (i.e., the lower and upper limits corresponding to 
Lm/Ls = 5 and 15 respectively). 

DOGs thus provide a significant contribution to the 
IRLD at z w 2. According to Caputi et al. (2007), 
the total IRLD from all 24/im-selcctcd z ~ 2 galax- 
ies is log(IRLD To tai) ~ 8.819±°;oyg, with the contribu- 
tion from ULIRGs alone being 42+22% °f tne total. Al- 
though DOGs presumably represent only a fraction of 
the z ~ 2 population (i.e., the reddest subset), they ap- 
pear to contribute s=s 26±14% of the total IRLD from the 
24^m-selected z ~ 2 galaxy population, and w 60^% 
of the IRLDulirg- The Caputi et al. (2007) z ~ 2 
sample is identified (in the GOODS-N and -S fields) us- 
ing photometric redshift estimates and should, in prin- 
ciple, include DOGs. What is noteworthy, however, is 
that the DOGs (selected solely on the basis of extreme 
R — [24] color) comprise such a significant component of 
the IRLDulirg i which suggests that the DOG selection 
criterion effectively selects out the bulk of ULIRGs at 
z ~ 2. 

Not all z <~ 2 galaxies may be selected as 24/xm sources. 
Indeed, one of the most successful methods of select- 
ing high-redshift galaxy populations is based on the UV 
emission of young, star- forming galaxies (e.g., Reddy 
et al. 2007, and references therein). Recently, Reddy 
et al. (2007) have estimated the contribution of UV- 
selected galaxies to the IRLD at z ~ 2. The IRLD 
contributed by DOGs is 158%, 60%, 35% and 56% of the 
IRLD contributed by the L m = 10 9 " 10 , 10 10 ~ n , 10 11 " 12 , 
and > 10 12 Lq UV-bright populations respectively. 

Do the bright and faint DOGs contribute compara- 
bly? To investigate this, we divided the DOG samples 
at F24 fim = 0.6 mJy, resulting in 2,149 galaxies fainter 
than the cut with a median flux density of 0.37 mJy, 
and 452 galaxies brighter than the cut with a median 



flux density of 0.85mJy. The redshift distributions on ei- 
ther side of this cut are slightly different (z s=s 1.9/2.1 for 
the fainter/brighter galaxies with spectroscopically mea- 
sured redshifts), but there are only 22 redshifts below 
the cut, so the distribution is not well constrained. If we 
assume that the distributions are the same at all flux den- 
sities, we find the IRLDs for the bright and faint subsam- 
ples are log[IRLD DO G(F24 M m > 0.6)] = 7.78 ± 0.02 and 
log[IRLD DO G(0.3 < F 24AI m < 0.6)] = 7.91 ± 0.01 respec- 
tively. The faint and bright sources thus contribute com- 
parable amounts to the IRLD, with the increasing num- 
ber making up for the lower luminosity of the fainter 
source population. 

In actuality, there is a large range in optical-to-infrared 
color, and we have chosen the (R — [24]) > 14 crite- 
rion simply to isolate the higher redshift galaxies. Since 
the redshift distribution is broad, selecting objects at a 
bluer threshold would still result in high-redshift pop- 
ulation, but with a larger low-rcdshift contamination. 
Moreover, we recall that the number density of DOGs 
rises toward fainter 24/^m flux densities, suggesting that 
i*24 < 0.3mJy sources may contribute significantly to 
the IRLD as well. In addition, we have not corrected 
our IRLD estimates for the sample incompleteness at the 
faint 24/im limit of our data. For all these reasons, the es- 
timates of space density and luminosity density reported 
here for the DOGs are likely to be lower limits, corre- 
sponding to the contribution of only the most extreme 
(in F 2 4 flm and (R— [24]) color) subset of the population. 
It is therefore remarkable that even this relatively rare 
population can contribute such a significant fraction of 
the overall infrared luminosity density. 

4.3.1. Comparison to Other z ~ 2 Galaxy Populations 

The galaxies selected by virtue of their extreme mid- 
infrared to optical flux density ratios seem to be under- 
represented in most optically selected samples of high- 
redshift galaxies. Indeed, our 24/im flux density crite- 
rion for DOGs (i.e., i^^m > 0.3 mJy, imposed by the 
existing MIPS observations of the Bootes Field) selects 
a relatively rare population that has only a handful of 
members in most existing small-field, deep samples. In 
this subsection, we attempt to compare our sample with 
other known z w 2 galaxy populations, most of which 
have been selected using observations at optical or near- 
infrared (i.e., rest- frame UV or optical) wavelengths. 

As we mentioned in §3.3.3 (and shown in Figures 10 
and 11), the DOG population is systematically fainter 
than the UV-bright populations selected at z ~ 2 by, say, 
the "BX/BM" selection criteria of Reddy et al. (2007). 
This is not surprising, since we are selecting galaxies to 
have a very large mid-infrared to UV flux ratio. What is 
interesting is that our selection results in a sample that 
lies at high redshift, with no bright, low-redshift mem- 
bers. Moreover, it is important to note that the DOG 
population, while being faint at optical wavelengths, is 
actually bolometrically very luminous, and thereby con- 
tributes significantly to (and perhaps dominates) the 
top end of the bolometric luminosity function (i.e., at 
L IR >1O 12 L ). 

A comparison with the z w 2 — 3 "Distant Red Galaxy" 
(DRG) samples selected using near-infrared selection 
techniques (typically, J — K > 2.3; e.g., Franx et al. 
2003; van Dokkum et al. 2004) is more difficult, since 
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the Bootes NDWFS field lacks very deep near-infrared 
imaging data. However, the DRG samples appear to 
have median Rab ~ 25.9 with 25(75)%-ile values of the 
distribution at Rab ~ 25.1(26.7) (van Dokkum et al. 
2006). The median rest-frame luminosity of DRGs is 
M 2200A ~ ~ 18 - 3 AB ma S for K < 21 (Vega) samples 
selected from the Hubble Deep Field South and MUSYC 
Deep fields (G. Rudnick, personal communication; Ga- 
wiser et al. 2006). These magnitudes, magnitude ranges 
and rest- frame luminosities are comparable to those mea- 
sured for our (R — [24] ) > 14 population, suggesting that 
these populations may overlap. It is now understood that 
the DRG selection results in a mixed sample of passively 
evolving systems, star-forming galaxies and AGN (e.g., 
Kriek et al. 2006a, b). Since our selection criteria require 
galaxies to be bright at 24/xm, it is likely that our selec- 
tion does not include the passively evolving galaxies. 

Recently, much work has been done using the BzK se- 
lection technique pioneered by Daddi et al. (2004) to se- 
lect star-forming and passive galaxy populations at high 
rcdshift (i.e., z « 2). We examined the (R— 24) colors of 
the BzK selected star-forming galaxies in the GOODS- 
N field (Daddi et al. 2007b,a), and find that only 4 of 
the 187 BzK galaxies in GOODS-N satisfy our selection 
criteria. The vast majority of the BzK galaxies are less 
extreme in their (R — [24]) colors, and are fainter in their 
24/im flux densities. This may be due to the lower AGN 
contribution (by selection) to the mid-infrared luminosi- 
ties of the BzK star-forming population. 

Given the surface density of bright DOGs (i.e., with 
> 0.3mJy, we would have expected roughly 14 
galaxies in the GOODS-N field. Indeed, there are 21 
galaxies with these selection criteria in the GOODS-N 
field, 3 of which are brighter than F 2 4 Mm = 1.0 mJy (E. 
Daddi and M. Dickinson, personal communication). Fif- 
teen of these 21 galaxies are also selected by the BzK 
techniques, but all but 4 are rejected from the BzK 
'star-forming' galaxy sample because they are hard X- 
ray sources (E. Daddi, personal communication). The 
(R — [24]) > 14 and BzK selections therefore appear to 
be comparable, with the star-forming systems at bright 
24/im flux densities effectively selected by both tech- 
niques. The simple (R— [24]) selection results in a slightly 
larger sample with perhaps a slightly broader redshift 
distribution. It is worth noting that the BzK technique 
can be fairly expensive in telescope time, requiring deep 
imaging in three bands. The simple technique described 
here of using the i?-[24] color selects similar star- forming 
galaxy and AGN populations more economically (90 sec 
with Spitzer/MlPS and 6000 sec of i?-band per field). 

Recently, a similar population of dust- 
obscured galaxies was identified by Yan et al. 
(2007) using the following selection criteria: 
F 2 ^m > 0.9 mJy; (i/F v (2AiJ,m)/uF v (8fim) > 
3.16; {uF v {2Anm)/vF v {Q.7^m) > 10. Our selection 
criterion for DOGs corresponds to 

{vF v {2Aiim)/vF v {Q.7iJ,m) > 28.6, and thus selects the 
redder, and a higher rcdshift subset of the population. 
The Yan et al. (2007) galaxies with (R - [24]) > 14(15) 
all lie at z > 1(1.6). The (^(24/xm)/^(8/zm) 
distribution for DOGs shows that the selection criteria 
used by Yan et al. (2007) would reject 36% of the DOGs. 
Finally, although dust-obscured galaxies have been 



found in other studies using comparable selection cri- 
teria, we emphasize the simplicity and robustness of the 
current approach. In particular, the simple selection cri- 
teria will make it easier, in principle, to quantify the 
selection function once a large sample of spectroscopic 
rcdshifts become available. One significant difference 
in the way we select the sample of DOGs is that we 
do not discriminate between AGN-dominatcd and star- 
formation-dominated systems. The mid-infrared colors 
of many of the DOGs suggest that they may indeed be 
dominated by AGN at bright 24/zm flux densities (i.e., 
^24^im <^ 0.8). Although the majority of the fainter pop- 
ulation does show stellar "bumps", PAH emission fea- 
tures in their mid-infrared spectra, and other evidence 
of starlight in their spectral energy distributions, their 
mid-infrared, far-infrared, or bolometric luminosity may 
yet contain a significant contribution from an AGN. The 
combination of AGN and star-formation signatures in the 
DOG population (and in some cases in the very same ob- 
jects) suggests that the population may be undergoing 
both rapid black-hole growth and rapid star-formation. 

4.4. The Possible Role of DOGs in Galaxy Evolution 

We have demonstrated that the simple selection cri- 
terion of large mid-infrared to optical flux density ratio 
efficiently isolates a significant population of z w 2 galax- 
ies. These systems are extremely luminous at rest-frame 
mid-infrared wavelengths, presumably due to warm dust 
heated by some combination of AGN and hot, young 
stars. What is the role of these galaxies in galaxy evolu- 
tion? 

In a seminal paper, Sanders et al. (1988) proposed that 
local ULIRGs represented a brief dust-enshrouded stage 
in the formation of quasars. Despite the extreme colors of 
the DOGs and the lack of local galaxies (even ULIRGs) 
with similar colors, there are some similarities between 
the two populations. They are both at the extreme lu- 
minous end of the galaxy population at their epochs, 
and they both appear to be heavily enshrouded systems 
exhibiting evidence for both intense star formation and 
AGN accretion activity. 

We therefore speculate that the DOGs represent a 
fairly short-lived phase in the evolution of massive (i.e., 
<~ AL*) galaxies, in between the phases represented by 
the sub-millimeter galaxy population (SMGs) and the 
more "passive" galaxies (e.g., selected as passive BzKs or 
DRGs). We envision a scenario in which gas accumulates 
in deep potentials, triggering an early episode of star- 
formation which quickly results in the formation of large 
quantities of dust. This results in a system which is lumi- 
nous at sub- mm wavelengths (as an SMG), with cold dust 
temperatures. The mid-infrared spectra of these systems 
would show PAH features, characteristic of star-forming 
systems, as has been observed (e.g., Menendez-Delmestre 
et al. 2007; Pope et al. 2008). At some point, as the ac- 
cretion and star-formation proceed, an AGN is triggered 
which heats the dust to warmer temperatures. It is at 
this point that the system would have a warmer char- 
acteristic dust temperature and be selected as a DOG. 
If the AGN is able to destroy or expel most of the dust 
(or the star-formation proceeds to the point where the 
dust and gas are sufficiently consumed), the DOG would 
evolve into an optically visible AGN. Whether or not the 
DOG is transformed into an optically-luminous AGN will 
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depend on the relative timescales of dust destruction / 
consumption and AGN accretion. It has been suggested 
recently that AGN feedback and star-burst driven winds 
can provide two important mechanisms for terminating 
the star-formation in galaxies (e.g., Hopkins et al. 2007, 
and references therein). In such a scenario, the galaxy 
evolutionary phase represented by DOGs would follow 
the bulge growth phase (i.e., the SMGs?), but precede 
the phase when the galaxy may be visible as an optical 
quasar or, eventually, as a red, passively-evolving galaxy 
(i.e., the DRGs or passive BzK galaxies?). If, instead, 
the AGN phase is more short lived than the UV-bright 
star-forming phase, then DOGs would evolve into a UV- 
bright star-forming population. 

There is some circumstantial support for such a sce- 
nario from the measured space density and redshift dis- 
tribution of the DOGs. The DOG surface density on the 
sky (of 0.089 arcmin~ 2 ) is similar to that of the lumi- 
nous sub-mm galaxy population with 850^m flux den- 
sities F 85 o Mm > 6mJy (Coppin et al. 2006), and they 
have a comparable redshift distribution (Chapman et al. 

2005) . The 850/xm sub-mm flux densities of the DOG 
population have not yet been measured, but the upcom- 
ing SCUBA-II Legacy survey should provide useful con- 
straints. In the scenario described above, the DOGs are 
expected to have warmer dust temperatures than the 
SMGs, and it is indeed interesting that the detection 
rate of SMGs at MIPS wavelengths (24^m and 70^m) 
is low, with only <40% of the SMGs in the GOODS-N 
field having 24/xm flux densities F24 flm > 0.3mJy (Pope 
et al. 2008). The space density of the brightest DOGs 
(with -F24/im > 1 mjy) is comparable to that of unob- 
scured QSOs at comparable redshifts (e.g., Brown et al. 

2006) . If these brighter sources are powered primarily 
by AGN, they represent a significant component of the 
accretion history of supermassive black holes, missed by 
traditional QSO surveys. 

Clustering analyses of the population may provide a 
clue to the possible evolutionary products of DOGs and 
the lifetime of the DOG phase relative to the SMG phase, 
and we will investigate this in a future paper. 

5. SUMMARY 

KPNO i?-band and Spitzer MIPS 24/im surveys of the 
Bootes Field of the NOAO Deep Wide-Field Survey re- 
veal a population of galaxies with very extreme opti- 
cal to mid-infrared colors. Galaxies with such red col- 
ors are absent from the local universe; e.g., there are 
no galaxies with these extreme (rest-frame) flux density 
ratios among the z < 0.1 galaxies in the Bootes Field. 
From the 8.14 deg 2 24/zm/i?-band overlap area, we se- 
lect a sample of 2,603 galaxies with 24 /jm flux densi- 
ties F 2 4 ftm > 0.3 mJy and (R - [24]) > 14 Vega mag 
(corresponding to i*24, tm /i<o.7fim ^ 1000. These galax- 
ies increase in number and as a fraction of the total 
number towards fainter 24/um flux densities, rising from 
7±0.6% of the population at i^m > ImJy to 13±1% 
at F24 flla w 0.3mJy. This reddening of the overall source 
population and the rise of the red fraction is not due to 
a simple fc-correction; the colors of this population of ob- 
jects are redder than any local ULIRGs, suggesting that 
these are very dust enshrouded objects which constitute 
an increasing fraction of the source population to fainter 
24^m flux densities. 



Selection of the reddest sources using a simple color 
cut of (R — [24]) > 14 effectively isolates a high red- 
shift population. The spectral energy distributions of 
these objects range from power-law SEDs suggestive of 
highly obscured AGN to SEDs exhibiting a "bump" at 
mid- infrared wavelengths. The fraction of sources with 
"bump" SEDs increases toward fainter 24/im flux densi- 
ties. We interpret the "bump" in the SEDs of the latter 
group as arising from a stellar component peaking at a 
rest- frame wavelength of 1.6/zm. 

Using spectrographs on the Spitzer Space Telescope 
and the telescopes of W. M. Keck Observatory, we have 
measured spectroscopic redshifts for 86 galaxies and have 
determined that the majority of this population lies be- 
tween the redshift range 1.5<z<2.5. The redshift dis- 
tribution may be modeled as a Gaussian with a mean 
redshift of 1.99±0.05 and width a = 0.45 ± 0.05. The 
high redshifts imply prodigious luminosities for this pop- 
ulation. Corrections to bolometric luminosities are un- 
certain, but crude estimates suggest luminosites of <~ 
10 12 ~ 14 L Q . Although the dominant source of the lu- 
minosity is not conclusively established, it is very likely 
that both AGN and prodigious star-formation contribute 
to the luminosity. These dust-obscured galaxies show a 
range in mid-infrared color, probably due to a change 
in the relative dominance of AGN over star formation 
contributions to the mid-infrared luminosity. A rough 
estimate of the space density of this population suggests 
a comoving space density of2.82x 10~ 5 /i 3 Mpc~ 3 . This 
dust obscured galaxy population contributes a signifi- 
cant fraction of the infrared luminosity density at z w 2, 
constituting roughly log IRLD sa 8.228 + 0.176 - 0.302, 
where IRLD is in units of L Q Mpc~ 3 . This is roughly 
26±14% of the total IRLD contributed by galaxies at 
this redshift, and 601^°% of the IRLD contributed by 
luminous infrared galaxies with Lir > 1O 12 L0. 

The properties of this population suggest that they 
may be galaxies caught in the very process of formation, 
undergoing rapid growth of their stellar components and 
/ or nuclear black holes. We speculate that these dust- 
obscured galaxies (or DOGs) may represent an evolu- 
tionary phase in the formation of massive galaxies, when 
the AGN is turned on and begins the process of ter- 
minating star-formation. Such a picture would envision 
DOGs as being the intermediate stage between the sub- 
mm galaxy population and the more quiescent galaxies 
(e.g., the 'passive' DRGs). 
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